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and in fact become more like the flux difference splitting. Since the differentiability and 
simplicity are desirable properties, one would still wish to search for a better splitting 
scheme that is strictly based on the flux vector splitting. 

In this paper, we propose a family of higher-order polynomial expansions for the mass 
flux that diminishes the diffusion error as M 0. We give a detailed study of the 
accuracy of the scheme for 1-D conical flow and 2-D shock wave/boundary-layer inter- 
actions. The weakness of the scheme is also pointed out and possibile improvements 
suggested. 


Analysis 

To exemplify the concept, let us consider the quasi two-dimensional system of equations 
for conical flows: 

dV dF „ 
dt + 8 rj ~ 

where U T = (p,ptx, pv y pE), F T = (pVypvu y pv* + p,pvH) } E = e + l/2(u 3 + v 2 ), and 
H = E + p/p. The flow considered consists of a very thin shear layer at the wall 
and a shock wave away from the wall. An algorithm must be capable of minimizing 
the numerical smearing(diffusion) at the locations where an eigenvalue changes sign or 
approaches zero. For example, Van Leer’s splitting[l] can represent shock profile well, 
while greatly diffusing the boundary layer. The Van Leer split mass fluxes are: 

Fi = F* + F ~ ; Ft = ±pa/A(M ± l) 2 . 

The net difference from the curve it approximates is largest at M = 0; its value equals 
pa/2. This error, viz numerical diffusion, significantly broadens the boundary layer, 
leading to incorrect velocity and temperature profiles. A simple way to remove the 
diffusion at M = 0 is by adding an extra higher-order term that allows the split mass 
fluxes to pass through the origin(Fig. 1), i.e., 

Ft = ±po/4[(M ± l) 2 + mi(M)(M 3 - l) 2 ], 

where the higher-order term has a coefficient mi, in general function of M. It should 
have the following properties: 

(1) n»i — ► —1 as M — * 0; 

(2) mi(Af) = 

(3) mi — * 0 as M — * ±1. 

A formula satisfying those properties is chosen as: 

m 1 = (Af 3 — l)/(Af 3 + l) s , 

where the exponent 5 is a free parameter; also shown in Fig. 1 is m t vs M with 5 = 4. 
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Fig. 1. Mass Flux Splitting in HOPE. Fig. 2. Pressure Splitting. 

In the conical flow calculations, the accuracy and convergence appear to be insensitive 
to the specified values of S = 2,4, 6. Now, regarding the flux aB a sum of convective and 
pressure terms, we can write the splitting formula for the flux vector: 




= F? 




With the realization in [5] that the pressure splitting could be considered separately in 
t e Van Leer formula[l], a whole host of freedoms for the pressure splitting becomes 
possible. Following is the list of formulas tested: 


(pi) : p* = =pl/4 (M ± 1)\M ^ 2 )p, 

(p2) : p* = (pi) T 3/4 M(M 2 - 1 ) a p, 

(p3) : P* = (pi) ± 3/4mi M (Af 2 ~ l) a p, 

and 

(p4) : P* = 1/2(1 ± 7 M)p. 

Figure 2 displays the distribution of the split pressure vs Af . The first formula is that 
used by Van Leerfl]. The second and third splits, (p2) and (p3), yield vanishing 
slope at M = 0, thus corresponding to central differencing. However, no instability 
was encountered in the conical flow problem with the (p2) or (p3) split used in an 
implicit code. The fourth split (p4) is obtained from an approximate integration along 
characterics. As will be seen later, the four formulas give essentially the same results 
for the conical flow calculated. 

Results And Discussion 

In this paper, two cases were tested to check the accuracy and convergence of the 
HOPE scheme. The first case is the 1-D self-similar conical flow over a 10-degree 
half cone at hypersonic speed, for which a detailed comparison study was conducted. 
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The flow conditions are: = 7.95, and Re x = 4.2 x 10*. Since P r = 1.0, exact 

solution gives adiabatic wall temperature, 13.64TOC . The second case is the 2-D shock 
wave/laminar boundary-layer interactions, for which experimental measurements were 
available[6]. The conditions are: M„ = 2.0, Re M = 2.96 x 10‘, and oblique shock 
angle /? = 32.585 degrees. In both cases, the results from the Roe splitting are also 
included for comparison. An implicit Newton iteration procedure was used to achieve 
steady-state solution with residual dropped by five orders of magnitude. 

Figures 3 and 4 show the pressure and temperature distributions from the first- and 
second-order solution on a 65-grid; little difference is seen. A monotone solution across 
the shock is obtained with the first-order scheme while oscillation appears in the second- 
order scheme, which can be eliminated by a TVD procedure. It is noted that the first- 
order pressure is smooth at the edge of the boundary layer, unlike the Roe solution which 
shows a slight discontinuity(not shown here). Although the boundary layer exhibits a 
steep temperature gradient, the HOPE scheme predicts the wall temperature correctly, 
indicating removal of the numerical diffusion associated with the original Van Leer 
splitting. 
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Fig. 3. Pressure Profile of Conic Flow. Fig. 4. Temperature Profile of Conic Flow. 


Figure 5 displays the results using various pressure splittings; they are practically iden- 
tical except the Van Leer pressure split (pi) shows some minor oscillation near the wall. 
However, the pressure splittings show significant effect on the convergence rate. The 
(p3) and (p4) splits are the best, comparable to the Roe splitting, while the other two 
are roughly two to three times slower. These may indicate possible instability in a more 

complex case. . 

Finally, for the 2-D case, the surface pressure and friction coefficient are plotted in 
Figs. 7 and 8. The first-order HOPE results compare fairly with Roe’s splitting and 
experimental data. However, the second-order calculation experienced difl&culty in con- 
vergence in which the residual was reduced by only two orders of magnitude and the 
result is not presented here. 
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Fig. 5. Comparison of Pressure Splits. Fig. 6. Convergence History. 


We suspect that a further investigation on other pressure splittings may lead to success 
in stability and convergence. Nevertheless, a systematic study of the eigenvalues of 
the split fluxes and the complete discretized system will prove to be a useful endeavor. 
Above all, the present research suggests that there are still possibilities in flux-vector 
splitting after Van Leer’s appeared nearly 10 years ago. The possibilities may very well 
still lie in the mass-flux splitting. 
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Fig. 7. Pressure at the Wall. 



Fig. 8. Friction Coefficient at the Wall. 
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